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Abstract. The Earth is a planet with a dipolar magnetic field which is agitated by a magnetized plasma wind 
streaming from the Sun. The magnetic field shields the Earth's surface from penetrating high energy solar wind 
particles, as well as interstellar cosmic rays. The magnetic dipole has reversed sign some hundreds of times over 
the last 400 million years. These polarity reversals correspond to drastic breakdowns of the dynamo action. The 
question arises what the consequences for the Earth's atmosphere, climate, and, in particular, biosphere are. It 
is shown by kinematic estimates and three-dimensional plasma-neutral gas simulations that the solar wind can 
induce very fast a magnetic field in the previously completely unmagnetized Earth's ionosphere that is strong 
enough to protect Earth from cosmic radiations comparable to the case of an intact magnetic dynamo. 
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1. Introduction 

Paleomagnetic records show that the magnetism of Earth 
has reversed itself hundreds of times over the last 400 mil- 
lion years (Valet & Meynardier 1993; Juarez et al. 1998; 
Gee et al. 2000; Selkin & Tauxe 2000; Valet 2003). In 
fact, geomagnetic polarity reversals represent the most dy- 
namic feature of the Earth's magnetic field. The polarity 
reversals do not occur instantaneously. Rather, transition 
periods, that span some thousand years and are charac- 
terized by unstable varying magnetic fields with no clear 
shape, lay between the stable dipole field states. During 
the transition periods the magnetic field strength can drop 
well below 10% of the average value (Juarez et al. 1998; 
Guyodo & Valet 1999; Selkin & Tauxe 2000) which sig- 
nifies a drastic breakdown of the Earth's dynamo. In the 
present time the magnetic south pole has wandered over 
more than 1100 km during the last 200 years. The strength 
of the Earth's field has decreased by 5% per century. 
This decrease is by far the fastest that has been verified 
since the last total field reversal, the end of the so-called 
Matuyama chron, 730,000 years ago. Also, by statistical 
estimates the Earth's dynamo is overdue for a reversal. 
Thus, we have to expect a transition period character- 
ized by a very small Earth's magnetic field in the near 
future. Since the Earth's dipole field provides us with a 
shield against cosmic rays and solar high-energy radiation 
one may wonder about the consequences for life on Earth. 
Also, having in mind that Mars lost the atmosphere al- 



most completely after the final breakdown of the magnetic 



field |(Luhmann fc Bauer 1992) one may speculate that 
stripping by the solar wind could alter the Earth's atmo- 
sphere. Severe climate changes could result. Interesting 
enough, during the last Brunhes-Matuyama reversal, no 
major changes in plant and animal life have been de- 
tected. This may be partly due to the fact that the atmo- 
spheric layers block some fraction of the cosmic radiation 
by scattering. On the other hand, the cosmogenic radionu- 
clide production varies at least over the last 200,000 years 
as function of short-term variations of the magnetic field 



(Frank 2000) 



In this contribution, we consider the interaction of the 
solar wind with a completely unmagnetized Earth. When 
the solar wind encounters unmagnetized objects, such as 
Venus (Luhmann 1995) and comets (Konz et al. 2004)1 
magnetic barriers and ionopauses develop. Although the 
interaction of a fully ionized and a weakly ionized gas is 
very complex, an important characteristic can be iden- 
tified - the generation of magnetic fields caused by rela- 
tive plasma- neutral gas shear flows. It has been shown 
(Huba k Fedder 1993)| that this process operates in the 
Venus' ionosphere and is responsible for the non-dipole 
magnetic field measured there. The same process has been 
studied in detail for the generation of seed magnetic fields 
in emerging galaxies (Wiechen et al. 1998; Birk et al. 
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2002) and in circumstellar disks (Birk et al. 2003) A kine- 
matic estimate indicates that relatively strong magnetic 
fields are generated in the Earth's ionosphere. So far, 
the only way to study the dynamical non-linear interac- 
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tion of the magnetized fully ionized solar wind plasma 
and the partially ionized Earth's ionosphere are three- 
dimensional plasma-neutral gas simulations. Our numeri- 
cal studies show the draping of the magnetic field of the 
solar wind and the self-generation of a relatively strong 
magnetic field in the Earth's ionosphere. 

2. A kinematic estimate of the ionospheric 
magnetic field induced by the solar wind 

The interaction of a fully ionized plasma with a partially 
ionized gas can be described by the fluid balance equa- 
tions for the mass densities, momentum densities and ther- 
mal pressures of the different species (see Sec. 3) together 
with the generalized Ohm's law. Ohm's law connects the 
electric fields and electric currents in a plasma. For the 
low-frequency dynamics we are interested in, it can be de- 
duced from the inertialess electron momentum equation 
(Mitchner & Kruger 1973; Wiechen et al. 1998) 

= -Vpo - n c e + ^v c x B^j - p c ^ oi (v c - v ; ) - 

Pencil (v - v n ) (1) 

where E and B arc the electric and magnetic fields, respec- 
tively, and by p c , n c , p c , v a , and v a p the electron pressure, 
particle density mass density, bulk velocities of the differ- 
ent plasma components (electrons, ions, and neutrals) and 
collision frequencies are denoted. The constants e and c are 
the elementary charge and speed of light. The electron ve- 
locity v c can be written in the form v c = v — ra\ j/ep 
with the plasma mass density p = p c + pi w p u the plasma 
bulk velocity v = {p c v c + p^vC)/ p w v ; and the electric 
current density j = ne(vi — v ) (n = n c = rij is the par- 
ticle density of the quasi- neutral plasma). Consequently, 
the inertialess equation of motion of the electrons yields 
Ohm's law in the form 
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where Ampere's law has been used to eliminate the current 
density. By Faraday's law one obtains an equation that 
governs the dynamical evolution of the magnetic field 
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with the magnetic diffusivity n = c 2 ^ e im e / 'Anne 2 . The 
final term, that results from the collisional momentum 
transfer between electrons and neutrals due to shear flow, 
is crucial for the generation of the magnetic fields un- 
der investigation. The collision frequency v cn is given by 
(Huba 19 98) v cn = n D <T(kT c /m c )^ 2 with a scattering 



cross section of a ~ 4 x 10~ 15 cm~ 2 and k, T e denoting the 
Boltzmann constant and the electron temperature, respec- 
tively. The main interaction between the solar wind and 
the Earth's ionosphere is expected to happen at a iono- 
spheric height of about 350 km where the ram pressure of 
the solar wind with an average velocity of vsw = 450 km/s 
is balanced by the ionospheric pressure. The collision fre- 
quency of solar wind electrons with ionospheric neutrals 
can be estimated to be v cn k, 200Hz. The solar wind con- 
sists of a fully ionized plasma that carries along a magnetic 
field of about B sw ~ 2 x 10~ 4 G (see, e.g., |(Priest 1987) 
for the solar wind parameters). The shear length scale L, 
that defines an estimate for the nabla operator in the final 
term in Eq. 3, is the least known parameter in our anal- 
ysis. In fact, we scale our numerical results to this (free) 
parameter. Shear lengths of L = 10 km seem to be rea- 
sonable. The ionospheric height gives an upper limit for L. 
The order of magnitude of the generated magnetic fields 
in the rest frame of the neutrals can roughly be estimated 
from the final term of Eq. (3) as 



m c c vsw 



= 5 x 10~ 4 G 



200 Hz 



vsw 



L 



10 km 



(4) 



450 km/s 

where r is the generation time scale. The time scale of 
the diffusion of the magnetic field measured by n is much 
larger than r. An upper limit for r cannot be given from 
a pure kinematic consideration. The saturation magnetic 
field results from self-consistent simulations of the entire 
dynamics presented in Sec. 3. 

For the parameters given, a field strength compara- 
ble to the present dipole value is generated after only ten 
minutes in the ionosphere. Thus, magnetic fields can be 
generated very efficiently around the unmagnetized Earth. 

3. The numerical model 

The interaction of the solar wind with the Earth's iono- 
sphere can be modeled by a plasma-neutral gas two fluid 
description. In our simulations, the following normalized 
plasma - neutral gas equations are numerically integrated 

|=-V-M (5) 

^ = -V-(v nPn ) (6) 

d 

— (pv) = -V-(/3vv)-Vp+(VxB)xB-p!/ pn (v-v n )(7) 
d 

^(^V„) = -V • (p n V„V n ) - Vp n - / O n ^np(v„ - v) (8) 

(9B 
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where p = p c + Pi and 7 denote the thermal plasma pres- 
sure and the ratio of specific heats (all quantities are now 
to be taken as dimensionless). The index n identifies the 
neutral gas quantities. The effective plasma-neutral gas 
collision frequency is u pn — (m c is cn + mj^i n )/(m c + r/i;) 
where v pn p = f np p n holds to guarantee momentum con- 
servation. In the induction equation (Eq. 9), P en denotes 
the normalized expression for i> cn cm c /e. The final term 
in the induction equation results in self-magnetization 
caused by relative sheared plasma-neutral gas flows (Lesch 
et al. 1989; Huba & Fedder 1993; Wiechen et al. 1998). 
Normalization of the physical quantities is carried out 
by the typical plasma mass density, the typical length 
scale, the plasma Alfven speed and time scale. The ther- 
mal pressure is normalized to the magnetic pressure. The 
thermal pressure is measured by a typical number den- 
sity and the kinetic temperature. The source terms in the 
balance equations for the mass densities, momenta and 
thermal pressures are treated in a flux-corrected manner 



(Birk & Otto 1996) For the simulation shown, we use a 
numerical domain of 7 Earth radii around the Earth in 
each direction (with the exception of 3 Earth radii in the 
solar direction) and one million grid points. The Earth is 
located at x — y = z = 0, the solar wind enters the nu- 
merical box from the upper y-boundary. At all the other 
boundaries outflow conditions are chosen. The neutral gas 
of the ionospheric layer around the Earth shows a cosh~ 2 - 
decay toward the interplanetary space. 

4. Numerical Results: Draping and 
self-magnetization 

Fig. 1 shows an arrow plot of the solar wind velocity field 
at three different times (t=30 s; upper plot, t=60 s; middle 
plot and t=600 s; lower plot). The wind encounters the 
Earth and is deflected around the planet. 

The magnetic field lines carried by the solar wind are 
draped around the Earth (Fig. 2). The draping leads to an 
amplification of the magnetic field near the Earth by one 
order of magnitude. This effect is well known, e.g., from 
investigations on the interaction of the solar wind with the 
unmagnetized Venus (Russel 1993; Luhmann 1995). 

Close to the Earth, the momentum transfer between 
the charged particles of the solar wind and the neutrals 
of the Earth's ionosphere becomes important (see final 
terms in Eqs.|2||SJ). Consequently, a new strong non-dipole 
magnetic field is generated by the sheared relative plasma- 
neutral gas motion (see final term in Eq. EJ) . Parameter 
studies show the kinematic finding (see Sect. 2) that the 
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Fig. 1. The solar wind flow around the Earth in the z = 0- 
plane at t=30 s (upper plot), t=60 s (middle plot) and 
t=600 s (lower plot). 



strength of the generated magnetic field depends on the 
shear length L. The maximum shear length is fixed in the 
simulation by an appropriate choice of the profile for i> cn . 
For a shear length L = 10 km a magnetic field of about 
the present dipole strength (S max ~ 0.3 G) is induced 
in the ionosphere after about 10 minutes(Fig. 3). For a 



magnetized Earth 




Fig. 2. The draping of the magnetic field lines of the solar 
wind around the Earth at t=600 s. 

given L the time scale of the field generation r results 
from the dynamics. The field is generated in heights of 
some hundreds of kilometers all around the Earth with 
the exception of the subsolar region where the magnetic 
field is weaker. If the shear length were chosen as say 100 
km, the maximum of the generated field strength would 
be B max « 0.03 G. 

We find that the draping effect is much weaker than 
the magnetic field self-generation by the shear flow. 

5. Conclusions 

We studied the interaction of the magnetized fully ion- 
ized solar wind plasma with the unmagnetized partially 
ionized Earth's ionosphere. When the solar wind hits the 
Earth the magnetic field lines carried along with it are 
draped around the planet. What is more important, the 
relative motion between the solar wind plasma and the 
ionosphere results in the self-generation of magnetic fields 
in the ionospheric layer. The strengths of these fields de- 
pend on the shear length of the relative flows, which, in 
contrast to the other relevant physical parameters, is not 
well known. For a reasonable shear length of L = 10 km 
the maximum strength of the newly generated magnetic 
field is comparable to the one of the present dipole field. 
Consequently, even in the case of a complete breakdown of 
the Earth's dynamo, the biosphere is still shielded against 
cosmic rays, in particular coming from the sun, by the 
magnetic field induced by the solar wind. 




Fig. 3. The strength of the induced magnetic field in the 
Earth's ionosphere at t=500 s. 
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